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ABSTRACT 

Spiral shocks are potentially a major source of turbulence in the interstellar medium. To address this problem 
quantitatively, we use numerical simulations to investigate gas flow across spiral arms in vertically stratified, 
self-gravitating, magnetized models of galactic disks. Our models are isothermal, quasi-axisymmetric, and 
local in the quasi-radial direction while global in the vertical direction. We find that a stellar spiral potential 
perturbation promptly induces a spiral shock in the gas flow. For vertically stratified gas disks, the shock 
front in the radial-vertical plane is in general curved, and never achieves a steady state. This behavior is in 
sharp contrast to spiral shocks in two-dimensional (thin) disks, which are generally stationary. The non-steady 
motions in our models include large-amplitude quasi-radial flapping of the shock front. This flapping feeds 
random gas motions on the scale of the vertical disk thickness, which then cascades to smaller scales. The 
induced gas velocity dispersion in quasi-steady state exceeds the sonic value for a range of shock strengths, 
suggesting that spiral shocks are indeed an important generator of turbulence in disk galaxies. 
Subject headings: galaxiesTSM — ISM:kinematics and dynamics — turbulence — instabilities — MHD 
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1. INTRODUCTION 

Turbulence in the interstellar medium (ISM) is observed to 
be pervasive and highly supersonic. Shocks created by ran- 
dom gas motions produce a rich variety of structures in the 
diffuse, atomic ISM as well as in gravitationally-bound giant 
molecular clouds (GMCs). Recent work has shown that turbu- 
lence is crucial to control of star formation within GMCs, and 
it is also believed to affect GMC formation processes. Tur- 
bulence thu s regulates star formation on both local and global 
scales (e.g.. iBallesteros-Paredes et aL "200^. In the absence 
of driving, shock dissipation and nonlinear cascades cause 
turbulence to decay on a time scale comparable to flow cross- 
ing times even for a medium with equipartition-strengt h mag- 
netic fields (e.g.. Stone, Ostriker, & Gammie 1998; M ac Low! 
ll999tlPadoan & Nordlund 1999), amounting to a few tens of 
Myr for the diffuse ISM. This implies that the ISM must be 
continuously stirred by one or (likely) more driving sources. 

A number of mechanisms have been proposed to drive tur- 
bulence, including H II region expansion, supernova explo- 
sions, and fluid instabilities involving magnetic fields and 
gravit y (see Mac Low & Klessen 2004; Elmegreen & Scalo 
12904' for recent reviews). Although less wefl recognized, 
galactic spiral shocks are also a n app ealing means to gener- 
ate ISM turbulence. Woodward ( fT976l) showed that cloud de- 
formation by passage of a shock, and the associated Kelvin- 
Helmholtz instabiliti es, together trigger random motions in 
the cloud. Wada & Koda (2004) suggested that an in-plane 
"wiggle instability" of spiral shocks can also drive turbu- 
lent motions, although this appears t o be suppressed when 
three-dimensional effects are included ("Kim & Ostrikei*2006', 
hereafter Paper I). Very recently, Bonnell et al. (2006) and 
|Dobbs et al. (2006) demonstrated that passage of a distribu- 
tion of clouds through a spiral shock gives rise to internal tur- 
bulent motions that follow Larson's (1981) empirical scaling 
law fairly well. 

Most of the studies cited above investigate two-dimensional 
dynamics driven by spiral shocks, neglecting the vertical de- 



gree of freedom. When the vertical dimension is taken into 
account and well-resolved. Paper I showed that non-steady 
motions of the spiral shock develop.' Non-steady behav- 
ior of spiral shocks and generation of vertical motions has 
in fact been seen in p revious n umerical models by other au- 
fliors (e.g., Martos & Coxlll998 : Gome z & Coxl l2002. 2004 
iBolev &'Durisen>.200 ^. In this Letter, we clarify the phys- 
ical causes of non-stationarity for spiral shocks in stratified 
disks. We also quantify the level of the induced gas motions 
to show that vertical spiral shocks should indeed be an impor- 
tant source of turbulence in spiral galaxies. 

2. METHODS AND MODEL PARAMETERS 

We study evolution of vertically-stratified gas flow across 
spiral shocks in local regions of self-gravitating, differentially 
rotating, magnetized galactic disks. Our studies use a mod- 
ified version of the ZEUS code (S tone & Normanl[T99 2a'bh 
to solve the time-dependent magnetohydrodynamic (MHD) 
equations presented in Paper I. We use the same local for- 
mulation described in Paper I, which explored fully three- 
dimensional models. Here, we perform much higher reso- 
lution models to focus on the detailed turbulent response, but 
suppress the degree of freedom parallel to the spiral arm (i.e. 
our models are quasi-axisymmetric). The reader is referred 
to Paper I for a complete description of our model prescrip- 
tion and numerical methods. Here, we briefly summarize the 
coordinate system and model parameters we adopt. 

We consider a local region centered on a tightly-wound 
spiral arm (pitch angle / <C 1), with a potential perturbation 
due to the stars assumed to be rigidly rotating with pattern 
speed ftp. We introduce a local, co-rotating Cartesian frame 
centered on {R,<j),z) = iRo,^pt,0)- The local frame is tilted 
such that X and y denote the in-plane directions perpendicular 
and parallel, to the arm, respectively, while z is perpendic- 
ular to the galactic plane (Roberts 1969). We set up a two- 
dimensional simulation domain with size x in the x-z 
plane (hereafter XZ plane), and assume all physical variables 
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TABLE 1 

Summary of model parameters and simulation results 



Model 


Go 


/3o 


F 


Grid 


ff(pc) 








Meff 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


~(9) 


(10) 


A 


1.8 


oo 


5 


1024^ 


196 


0.66 


0.63 


0.38 


4.0 


B 


1.8 


oo 


5 


512^ 


196 


0.67 


0.62 


0.36 


3.9 


C 


2.0 


oo 


7 


512^ 


218 


1.01 


0.86 


0.45 


5.0 


D 


2.5 


oo 


10 


5122 


272 


1.55 


1.20 


0.58 


6.2 


E 


1.5 


10 


5 


512^ 


169 


0.64 


0.52 


0.31 


3.8 


F 


1.8 


10 


7 


512- 


203 


0.93 


0.79 


0.42 


4.4 


G 


2.0 


10 


10 


5122 


225 


1.54 


1.19 


0.52 


5.4 



Note. — Col. (l)-(4): Model name and input parameters. Col. (5): Numerical resolution. Col. (6): Vertical scale height of gas in initial magnetohydrostatic 
equilibrium. Cols. (7)-(9): Density-weighted late-time velocity dispersions. Col. (10): Effective Mach number of the large-scale spiral shock. 



to be independent of y. The simulation box size is equal to 
the arm-to-arm separation. We allow for nonzero values of ve- 
locity V, and magnetic field B,.; initially v,, = {Q,Q-Vtp)RQ-VlQX 
and By is independent of x, where fio is the local angular ve- 
locity at /?o. We also include Coriolis forces. This quasi- 
axisymmetric approximation prevents clump-forming, non- 
axisymmetric gravitational instabilities from developing, en- 
abling us to focus on the properties of turbulence independent 
of cloud formation. 

We idealize the ISM by treating it as an isothermal gas with 
effective speed of sound c, = 7kms"'. The initial gaseous 
disk (without the spiral potential) is in vertical magnetohydro- 
static equilibrium. As reference values, we take Rq = lOkpc, 
fio = 26kms~' kpc~', and epicyclic frequency kq = 2'/^51() 
(for a flat rotation curve). The corresponding orbital period 
is forb = 27r/fi() = 2.4 X lO'^yr. For the spiral arm parameters, 
we adopt ftp = Q,q/2 and sin / = 0. 1 . We adopt = ttRq sin / = 
3.14 kpc for a two-armed spiral. We impose reflection sym- 
metry with respect to the midplane, and apply open boundary 
conditions at z = = 4H, where H (~ 170-270 pc) is the 
initial disk scale height. We adopt sheared-periodic boundary 
conditions in x. 

The three key parameters that characterize our model disks 
are 

i^QCs cl 2 /^|$sp| 



sin; 



(1) 



ttGSo 

where Eo is the gas surface density, va is the Alfven speed, 
and $sp is the amplitude of the imposed sinusoidal potential 
(Paper I). The quantity F measures the ratio of the perturbed 
sinusoidal radial force to the mean axisymmetric gravitational 
force (Roberts 1969). We do not consider the vertical varia- 
tion of the spiral potential perturbation; the perturbed vertical 
force is negligible (< 4%) compared to the background verti- 
cal force. In our simulations, the potential amplitude slowly 
increases from to F over ^ l.Sforb- We present results for 
seven numerical models, as listed in Tabled 

3. NONLINEAR SIMULATIONS 

3.1. Spiral Shock Evolution in a Stratified Disk 

In this subsection, we describe the evolution of spiral 
shocks in our fiducial model G; evolution in other models is 
similar. As the amplitude of the stellar spiral potential grows, 
magnetosonic waves emerge and steepen, forming a shock 
front near x/L^ = 0.04 by f/foib ~ 1-0. At this time, fluid 
motions are nearly horizontal, and the shock front is nearly 
vertical. The gravity due to the large post-shock midplane 
density enhancement soon breaks vertical force balance, how- 
ever, pulling material at high |z| toward the midplane. With in- 
creased postshock density, the shock near z = tends to move 



upstream toward the spiral potential minimum, causing the 
shock front to bend. Gas entering the curved shock is initially 
bent upward toward the shock front and then pulled downward 
due to strong vertical gravity. The falling material establishes 
repulsive pressure gradients and bounces back to high |z|. The 
instantaneous streamlines shown for model G at f / foib = 1 .2 in 
Figure[0? reflect gas motions at this stage. 

Our model simulations show that spiral shocks in the XZ 
plane never achieve a steady state, as illustrated by the strong 
density variations in Figure [Qj.c for model G. This result is 
quite unlike the steady solutions that obtain for models which 
neglect vertical degrees of freedom (e.g., IWoodwardlll975t 
Kim & Ostriker 2002). In particular, when the fluid vari- 
ables are allowed to vary with z, the gas leaving the right x- 
boundary does not in general have the same z as it had when 
it originally entered the left x-boundary^, as Figure [iji illus- 
trates. When the vertical gravity exceeds the vertical repulsive 
pressure force downstream from the shock, post-shock mate- 
rial is drawn toward the midplane, and the shock front at high 
|z| bends downstream (Fig.fy?). As gas is further compressed, 
vertical pressure gradients overwhelm vertical gravity and the 
gas even at \z\/H ~ 0.5 is able to rebound to \z\/H ^ 3. This 
in turn forces the shock front at large |z| to shift back upstream 
(Fig. [Q;). The vertically-rebounding gas overshoots equilib- 
rium, vertical gravity again dominates, and the cycle repeats. 
As gas traverses the spiral arms, therefore, the shocks in our 
models continue to quasi-periodically shift back and forth per- 
pendicular to the arms. These "flapping" motions of spiral 
shocks are strongest at high Random motions at scales 
comparable to the height H are also driven, and many weak 
shocks develop especially at high |z| (see Fig.^). 

In addition to large-scale flapping, we find that spiral shocks 
in model G experience two small-scale instabilities that aid in 
transfer of random gaseous kinetic energy to smaller scales: 
an advective-acoustic cycle and vortex generation. The 
advective-acoustic cycle occurs as the most overdense mid- 
plane region perturbs the upstream flows by launching acous- 
tic waves, analogous to the vortical-acoustic cycle for insta- 
bilities in isothermal Bondi-Hoyle-Lyttleton accretion (e.g., 
|R)glizzo 2002). Because the spiral shock front is curved in 
the XZ plane, Crocco's theorem ensures generation of vortic- 
ity. Moreover, strong rising and falling currents of the gas 
in the post-shock region create a sheared velocity field favor- 
able to Kelvin-Helmholtz instability, which can also gener- 
ate vorticity. We remark that vortices formed at the shock 
front are not fully resolved in our models because instabilities 

2 This is because the fundamental vertical oscillation period need not equal 
the arm-to-arm flow period, although these timescales are typically within a 
factor of a few of each other 
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Fig. 1. — Snapshots of density in logaritlimic color scale of model G at 
'/'orb = 1-2, 1.7, 1.8 from top to bottom. Solid lines in (a) represent instan- 
taneous streamlines of gas. The rectangular box A in (c) indicates a sector 
enlarged in Figurelll 

grow faster at smaller scales. Nevertheless, the velocity dis- 
persions are essentially independent of numerical resolution 
(see below), suggesting that these instabilities may enhance 
the small-scale cascade but do not affect the overall level of 
induced turbulence. 

Figure |2]plots typical velocity and density structures inside 
a spiral arm when turbulence is saturated. Many weak shocks 
and well-resolved vortices of both signs are apparent. The 
instantaneous, density-weighted, total velocity dispersion in 
the region shown in Figure |2] is 1.51 times the sound speed. 
We find that the power spectra of the perturbed velocities ex- 
tend smoothly over all scales, indicating fully-developed tur- 
bulence. Because of the large-scale shock in Vj, the large 
shear/streaming in v,,, and the vertical stratification, the power 
spectrum of v- as a function of best characterizes the true 
turbulence. At heights \z\/H < 1, this ranges from v?(fcv) cx 
-0.5 „ - LAx/ilir) < 8 to x^Ak,) oc kf -^ at 8 < n, < 16 

: modes. For 



/tj" " at n.v 
to v^(k^) oc k~ 

1 < \z\/H < 2, the power spectra are slightly steeper, having 



-1.3 



at 16 < < 64 for the k- 



v^fe) oc at < 8, vf{k,) oc k'' ' at 8 < n,- < 16, and 
vj{kx) oc A;"'-^ at 16 < n, < 64. 



3.2. Level of ISM Turbulence 

We now quantify the level of random gas motions driven 
by spiral shocks. Since the velocities in XZ spiral shocks are 
non-uniform and non-stationary, it is useful to construct the 
mean velocity field (v,) (with / = x, y, or z), where the bracket 
( ) denotes a time average over f/foi-b 4-8 after turbulence 
saturates. We then measure density-weighted velocity disper- 
sions using (jj = J pSvjdxdz/ J pdxdz, where (5v,- = v,- (v,).-' 
Figure|3jfl) plots cr,(r) for model G, while columns (7)-(9) in 
TablelUlist {afy^^ for all the models. 

' In computing each cr, , if we initially subtract out velocities corresponding 
either to unperturbed rotation or to a thin-disk spiral shock (instead of (v,)), 
then the velocity dispersions would be larger. 
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Fig. 2. — Velocity vectors and density in logarithmic scale at t/t„^i, = 3.2 in 
the region A marked in FigureQJ. The color scale is the same as in FigurelTI 
The size of the arrow above the box corresponds to the sound speed. 



As Figure ^a) shows, random gas motions in model G 
are supersonic in the x- and y-directions, and exhibit large- 
amplitude temporal fluctuations. The characteristic periods 
of these quasi-periodic fluctuations are in the range 0.5 - 
0.9forb- The maximum velocity dispersions occur when the 
spiral shock is temporarily vertical. Evidently, the system 
reaches a quasi-steady state in which dissipation of turbu- 
lence (in shocks and through cascades) is offset by the con- 
tinual input of new turbulent energy from the large-scale flap- 
ping. Turbulence is strongest in the immediate post-shock re- 
gion, and declines with increasing distance from the shock. 
Vertically-averaged velocity dispersions inside the spiral arms 
typically exceed those in interarm regions by about a factor of 
2. For the model parameters we have considered. Tabled in- 
dicates that (cTf) ((Ty) ^ 2(cr,). Comparison of (cr,) between 
models A and B shows that the velocity dispersion is indepen- 
dent of numerical resolution. 

What determines the level of velocity dispersions in XZ 
spiral shocks? Table [2 suggests that stronger shocks (hav- 
ing larger values of Qq\ and/or F) yield larger cr,. To 
characterize the shock strength in each model as simply as 
possible, we have run one-dimensional counterparts with the 
same set of parameters, using the thick-disk self-gravity pre- 
scription (see Paper I). We define the effective Mach number 
TWeff = {^i/^xY^^, where Si and E2 denote the preshock 
and postshock surface densities, respectively, in the result- 
ing (one-dimensional) spiral shock. Column (10) in Tabled 
lists A^eff for each model. Figure |3j^) plots the total veloc- 
ity dispersion CTtot = ((ctJ) + {(jy) + (cr?))'/^ as a function of 
yWeff- While cTtot monotonically increases roughly as CTtot/Ci = 
0.6A^eff- 1-4 for TWeff < 5.5, it is more or less constant at 
about 2cs for A^eff > 5.5. We have found that models with 
other values of fin/ilo follow the same crtot-A^eff relation as 
shown in Figure |3fe). 

4. DISCUSSION 

Galactic rotation can supply an effectively inexhaustible 
amount of kinetic energy to power turbulence in the ISM 
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Fig. 3. — (a) Time evolution of the density-weighted velocity dispersions 
(T, in model G. (h) Total velocity dispersion atot vs. the effective Mach num- 
ber A^eff of the spiral shock in models B to G. The errorbars represent the 
standard deviations in the temporal fluctuations of cr,. 

(Ifleck 198 1)."* In this Letter, we have shown that spiral shocks 
are in general non-stationary in the radial-vertical plane, and 
can efficiently transform some of the available bulk rotational 
energy into random gas motions. Vertical force imbalance 
drives radial flapping of the shock front at \z\/H ^ 0.2, which 
then feeds turbulence at smaller scales. The random gas mo- 
tions induced by the spiral shock persist despite strong shock 
dissipation, and yield time-averaged in-plane velocity disper- 
sions ^7-10 km s"' for a range of shock strengths, similar to 
the observed line widths of cold and warm atomic gas in the 
Milky Way (e.g.. iHeiles & Troland 2003). Vertical velocity 
dispersions are lower, but still amount to ~ 1 /2 of the thermal 
velocity dispersion. 

Although our presentation has focused on shocks with real- 
istic parameters, in fact turbulence generation by spiral shocks 
appears to require neither self-gravity nor magnetic fields. 
Our model simulations, including those unlisted in Table ^ 
show that anything (e.g., an imposed variation of $sp with 

To maintain a fixed level of turbulence, the energy dissipation per 
unit mass per unit time must be balanced by the net stresses, Q{pVy5vy - 
BxByliAiv)) I p (e.g. Piontek & Ostrikei 2005). If this power is supplied at 
the expense of overall ISM accretion toward the galactic center, the accre- 



z) that makes the primary shock non-vertical will end up pro- 
ducing a non-steady flow - simply because vertical oscillation 
periods need not (and in general do not) agree with the hori- 
zontal crossing time between arms. 

Since strong spiral arms usually imply a high rate of star 
formation, spiral shocks and stellar energy sources may work 
together in generating turbulence in many galaxies. The ab- 
sence of observed correlation between spiral arm phase and 
turbulent amplitude (e.g., Dickev et al. 1990), however, sug- 
gests that these two processes are not the only important 
sources of turbulence. In addition, radio observations of ex- 
tended H I disks in face-on galaxies show that the total vertical 
velocity dispersions are as large in the non-star-forming outer 
parts as in the star-forming inner regions (van Zee & BrvanJ 
again suggesting that additional sources of turbulence 
are able to compensate when needed. 

One compensating source of turbulence may be provided by 
the magnetorotational instability (MRI) ("Sellwood & Balbud 
|T999: Kim. Ostriker. & Stone 2003; Dziourkevitch et al] 
12004). The MRI should naturally be strongest exactly where 
spiral shocks are absent: in the outer galaxy, where gravity is 
weak, the disk flares, and the density drops; and in interarm 
regions, where the angular velocity decreases outward (in 
arms, shear is reversed so that MRI cannot occur). 

While in this Letter we adopt isothermal conditions for the 
gas and take the mean magnetic field parallel to the spiral 
arm, the real ISM has a multi-phase struc ture and is threaded 
also by weak vertical magnetic fields. iPiontek & Ostrikei 
(2004, 2005) demonstrated that MRI is able to generate ran- 
dom gas motions up to 8 km s"' in two-phase gas under fa- 
vorable conditions. It will be interesting to see how turbulent 
driving in spiral shocks develops for a multiphase medium, 
and whether MRI and spiral shocks can indeed establish a 
geographically-compensating balance of power. 
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tion time is ~ (2f2)"'p(f2S)^(pv,v(5vj)"', which exceeds the Hubble time if 

{pv,5vy-B,Byl{ATx))r^pc]. 
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